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' We propose how to realize a multiqubit tunable phase gate of one qubit simultaneously controlling 

n qubits with four-level quantum systems in a cavity or coupled to a resonator. Each of the n two- 
qubit controUed-phase (CP) gates involved in this multiqubit phase gate has a shared control qubit 
' but a different target qubit. In this propose, the two lowest levels of each system represent the 

two logical states of a qubit while the two higher-energy intermediate levels are used for the gate 
' (-H ' , implementation. The method presented here operates essentially by creating a single photon through 

^ the control qubit, which then induces a phase shift to the state of each target qubit. The phase 

shifts on each target qubit can be adjusted by changing the Rabi frequencies of the pulses applied 
, to the target qubit systems. The operation time for the gate implementation is independent of the 

. number of qubits, and neither adjustment of the qubit level spacings nor adjustment of the cavity 

^ ' mode frequency during the gate operation is required by this proposal. It is also noted that this 

approach can be applied to implement certain types of significant multiqubit phase gates, e.g., the 
multiqubit phase gate consisting of n two-qubit CP gates which are key elements in quantum Fourier 
transforms. A possible physical implementation of our approach is presented. Our proposal is quite 
general, and can be applied to physical systems such as various types of superconducting devices 
coupled to a resonator and trapped atoms in a cavity. 
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PACS numbers: 03.67.Lx, 42.50.Dv 

I. INTRODUCTION 

During the past decade, various physical systems have been considered for building up quantum information proces- 
sors [1-4]. Among them, cavity QED analogs with solid-state systems are particularly appealing [5-10]. Theoretically, 
it was predicted that the strong coupling limit, which is difficult to achieve with atoms in a microwave cavity, can 
be readily realized with superconducting charge qubits [11-13], superconducting flux qubits [14], or semiconducting 
quantum dots [15]. Moreover, the strong coupling cavity QED has been experimentally demonstrated with supercon- 
ducting qubits [16,17] and semiconductor quantum dots embedded in a microcavity [18]. These experimental results 
make solid-state qubit cavity QED a very attractive approach to quantum information processing. The goal of this 
work is to try to find ways to implement multiqubit gates in solid state QED qubit systems. Let us first introduce 
two-qubit gates and then generalize to multiqubit gates. 

A. Two-qubit controlled phase gates 

It has been shown that standard two-qubit controUed-phase (CP) gates, together with single-qubit gates, form the 
building blocks of quantum information processors. A standard two-qubit CP gate is described by the transformation 
|00) |00) , |01) |01) , |10) |10) , and |11) - |11). This shows that: when the control qubit (the first qubit) is 
in the state |1) , the gate fiips the phase of the state |1) of the target qubit (the second qubit) by tt (i.e., |1) — >■ — |1)) 
but does nothing otherwise. So far, a large number of theoretical proposals for realizing this two-qubit CP gate have 
been presented with many physical systems. Moreover, this two-qubit CP gate, together with a two-qubit controlled- 
not (CNOT) gate or a two-qubit «SWAP gate, has been experimentally demonstrated in, e.g., cavity QED [19,20], 
NMR [21], and superconducting qubits [22-25]. 

In this paper, we focus on a different type of two-qubit CP gate. To see the difference between this type of gate and 
the standard CP gate above, let us first consider a two-qubit CP gate with qubit 1 (the control qubit) and qubit k 
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FIG. 1: Schematic circuit of a multiqubit controlled-phase (CP) gate with one qubit (on top) simultaneously controlling n 
target qubits. The circuit on the left side is equivalent to the circuit on the right side. Each of the n two-qubit CP gates 
forming this multiqubit phase gate has a shared control qubit (i.e., qubit 1, represented by the filled circle located on top) but 
a different target qubit (i.e., either qubit 2, 3, ... , or n+ 1). The element 6k represents a controlled-phase shift by e'*'', which 
is performed on the target qubit k when the control qubit 1 is in the state |1) (fe = 2, 3, n+l). Namely, if the control qubit 
1 is in the state |1), then a phase shift e'**" happens to the state |1) of the target qubit k; otherwise nothing happens. 

(the target qubit). Here, k is the label of the target qubit, which is an integer greater than 1. The operator describing 
this gate is given by the following matrix 
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where < ^fe < 27r. Namely, if and only if the control qubit 1 is in the state |1), a phase shift e*^*" happens to the 

state |1) of the target qubit k. It is obvious that for 9k ^ n, this two-qubit tunable CP gate of Eq. (1) is different 
from the standard two-qubit CP gate. It is called tunable, because 9k can be adjusted, as described below. 

B. Multi-qubit phase gates 

Let us now consider a multiqubit quantum phase gate which consists of n two-qubit CP gates as depicted in Fig. 
1. From Fig. 1, one can see that the n two-qubit CP gates act on different qubit pairs (1, 2), (1, 3), and (1, n + l), 
respectively. Each two-qubit CP gate has a common control qubit (qubit 1) but a different target qubit (2,3, or 
n + 1). The CP gate acting on the qubit pair (l,fc) is described by the operator Ri^ above, which induces a phase 
shift e*^* to the state |1) of the target qubit k when the control qubit 1 is in the state |1) . In the general case, any two 



of the phase shifts , e 



, e''^"+i for the n two-qubit CP gates can be the same or different, and any one of them 



can be arbitrary. The operator characterizing this multiqubit phase gate with one qubit simultaneously controlling 

the n target qubits is 



where the operator Rik is given by Eq. (1). 

C. Motivation for studying multi-qubit gates 
There are several motivations for this work: 
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FIG. 2: (a) Schematic circuit of a multiqubit controllod-phasc gate with n control qubits acting on one target qubit located 
at the bottom of the circuit shown here. The clement Z represents a Pauli rotation CTz, i.e., a phase shift by tt around the 
^;-axis (with n control qubits on the filled circles). Namely, if the n control qubits are all in the state then the state |1) at 
Z is phcise-shifted by tt (i.e., |1) — > e^^ |1) = — otherwise nothing happens to the state at Z. (b) Schematic circuit of a 
multiqubit controUed-NOT gate with n control qubits acting on one target qubit. The symbol © represents a controUcd-NOT 
gate (with n control qubits on the filled circles). If the n control qubits are all in the state then the state at © is bit-flipped 
(i.e., |1) ^ |0) and |0) ^ |1)). 



(i) Quantum gates with multiple control qubits or multiple target qubits are of great importance in quantum 
information processing such as realizing quantum error-correction protocols, constructing quantum circuits, and im- 
plementing quantum algorithms. When using the conventional gate-decomposition protocols to construct a multiqubit 
controlled gate [26,27], the procedure usually becomes complicated (especially for a large n), as the number of single- 
qubit and two-qubit gates required for the gate implementation heavily depends on the number n of qubits, and 
therefore building a multiqubit gate may become very difficult. 

(ii) Recently, attention is shifting to the physical realization of multiqubit gates (e.g., [28]). Although several meth- 
ods for implementing multiqubit gates based on cavity QED or trapped ions have been proposed [29-32], realization 
of controllable multiqubit phase gates based on cavity QED or ion traps has not been thoroughly investigated. 

(iii) The methods previously proposed [29-32] are mainly for realizing quantum controUed-phase or controlled-not 
gates with multiple- control qubits acting on one target qubit (Fig. 2). These methods provide a simple and/or 
fast way for realizing quantum gates with multiple control qubits (a type of multiqubit gate significant in quantum 
information processing), when compared with the conventional gate-decomposition approaches. However, it is noted 
that these proposals cannot be extended to perform a different kind of important multiqubit phase gate, i.e., quantum 
phase gates with one qubit simultaneously controlling multiple target qubits. 

D. Advantages 

The purpose of this work is to find a way of implementing a multiqubit tunable phase gate of one qubit simultaneously 

controlling n qubits with four-level quantum systems in a cavity or coupled to a resonator. Our proposal operates 
essentially by creating a single photon through the control qubit, which then induces a phase shift to the state of each 
target qubit. As shown below, this proposal has the following advantages: 

(i) since the excited level |3) is unpopulated during the gate operation, decoherence due to spontaneous decay from 
this level is suppressed; 

(ii) neither adjustment of the level spacings of qubit systems nor adjustment of the cavity mode frequency during 
the gate operation is needed, thus the operation is much simplified; 

(iii) the operation time required for the gate realization is independent of the number of qubits and thus does not 
increase with the number of qubits; 

(iv) the phase shift on each target qubit can be adjusted by varying the Rabi frequencies of the microwave pulses 
applied to the target qubit systems; and 

(v) the n two-qubit CP gates forming this multiqubit phase gate can be simultaneously performed using six op- 
erational steps (independent of n). Thus, the gate operation is much simplified, especially when the number n is 
large. In addition, as shown below, the present method can be applied for implementing certain types of significant 
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FIG. 3: Illustration of four-level qubit systems. The energy eigenvalues for the four levels |0), |2), and |3) are denoted by 
Eo, El, E2, and E3, respectively. In (a), the level spacings satisfy E2 — Ei > Ei — Eo, E^ — E2; and E^ — E2 < Ei — Eq. In 
(b), the level spacings satisfy Ei — Eg > E2 — Ei > E3 — E2. In (c), the level spacings meet E2 — Ei > Ei — Eq, Eg — E2; 
and E'i — E2 > El — Eq. The four-level structure in (a) applies to charge-qubit systems, the one in (b) applies to phase-qubit 
systems, and the one in (c) applies to flux-qubit systems. In addition, the four-level structure in (b) is also available in atoms. 

multiqubit phase gates, e.g., the multiqubit phase gate consisting of n two-qubit successive CP gates, which are key 
elements in the quantum Fourier transform (QFT). 

We stress that this proposal is quite general, and can be applied to physical systems such as trapped atoms in 
a cavity and various types of superconducting qubit systems coupled to a resonator. We note that implementing a 
two-qubit tunable phase gate with two atoms based on cavity QED was previously proposed [33,34]. However, to the 
best of our knowledge, no one has yet demonstrated how to perform a multiqubit tunable phase gate with one qubit 
simultaneously controlling n qubits in cavity QED. Also, our work is the first to show how to perform the n successive 
two-qubit CP gates in QFT by the use of only one single-mode cavity. We believe that this proposal is useful since it 
provides a simple and general protocol for realizing a multiqubit phase gate for which the phase shift on each target 
qubit is tunable. 

This paper is organized as follows. In Sec. II, we briefly review the basic theory of four-level quantum systems 
coupled to a single- mode cavity and/or driven by classical pulses. In Sec. Ill, we show how to realize a multiqubit 
tunable phase gate with one qubit simultaneously controlling n target qubits, by the use of {n + 1) qubit systems in 
a cavity. In Sec. IV, we discuss how to apply the present method for implementing two types of significant quantum 
phase gates with multiple target qubits or multiple control qubits. In Sec. V, we give a brief discussion of the 
experimental feasibility for implementing a six-qubit controUcd-phase gate, which are key elements in producing a 
QFT. In Sec. VI, we compare our work with previous work. A concluding summary is presented in Sec. VII. 

II. SYSTEMS INTERACTING WITH A CAVITY AND PULSES 

The four- level quantum systems throughout this paper could be either natural atoms or artificial atoms (e.g., 
superconducting devices), which have the four levels shown in Fig. 3. Note that the four-level structure in Fig. 3(a) 
applies to superconducting charge- qubit systems [1], the one in Fig. 3(b) applies to phase-qubit systems [2,35], and 
the one in Fig. 3(c) applies to flux-qubit systems [1,36]. In addition, the four- level structure in Fig. 3(b) is also 
available in atoms. By adjusting the level spacings of the qubit systems or the cavity mode frequency, the level |0) 
can be made not to be affected by the cavity mode or the classical pulses (see the discussion given in the appendix). 
Under this condition, the level |0) can be dropped from the discussion below (also, see Fig. 4). We note that four 
levels are needed. This is because (as explained in Sec. Ill) a three-level quantum system is not enough to implement 
our desired gate. 

A. System-cavity-pulse resonant Raman coupling 

Let us consider a qubit system in a cavity or coupled to a resonator, which has a four-level structure depicted in Fig. 
3. The energy eigenvalues for the three levels |1), |2) , and |3) are denoted as Ei, E2, and E3, respectively. Assume 
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that the cavity mode is coupled to the |2) o |3) transition but decoupled (highly detuned) from the transition between 
any other two levels [Fig. 4(a)], which can be achieved by: (i) changing the cavity mode frequency, (ii) choosing the 
qubit system (e.g., atom) appropriately to have a desired level structure, and/or (iii) adjusting the level spacings of the 
qubit system. Note that the cavity mode frequency for both optical cavities and microwave cavities can be changed 
in various experiments (e.g., see, [37-41]). And, for superconducting qubit systems, the level spacings can be readily 
adjusted by varying the external parameters (e.g., the external magnetic flux and gate voltage for superconducting 
charge-qubit systems, the current bias or flux bias in the case of superconducting phasc-qubit systems and flux-qubit 
systems, see e.g. [1,2]). In addition, we assume that a classical pulse is applied to the qubit system, which is coupled 
to the [1) ^ |3) transition but decoupled (highly detuned) from the transition between any other two levels [Fig. 4(a)]. 
The Hamiltonian of the whole system can thus be written as 

3 

H = hiOca+a + ^Eill) hg{a~^a23 + il.c.) 

1=1 

+fil7(e"^*c7f3+H.c.), (3) 

where a+ and a arc the photon creation and annihilation operators of the cavity mode with frequency Uc', g is the 
coupling constant between the cavity mode and the )2) o |3) transition; is the Rabi frequency of the pulse, and w 
is the frequency of the pulse; ct^ = ]2) (3| , and = jl) (3| . 

Let us now assume that the cavity mode is off-resonance with the )2) <-> )3) transition, i.e., Aj, = lo^2 ~ ^ 9i 
and the pulse is off-resonance with the |1) O |3) transition, i.e., A = uj^i — w 3> f^, where 0732 = (-E3 — £^2)/^ 
and CJ31 = (i?3 — Ei)/h are respectively the |2) o |3) transition frequency and the |1) O |3) transition frequency 
[Fig. 4(a)]. Under these conditions, the level ]3) can be adiabatically eliminated [42]. Thus, for Ac = A^^, the 
effective Hamiltonian in the interaction picture is [43] 



Hi = -h 



^|l)(l| + |-a+al2)(2] + ^(a+<7++H.c.) 



(4) 



where crj*^ = |2) (1). The first two terms in Eq. (4) are ac-Stark shifts of the levels |1) and )2) induced by the pulse 
and the cavity mode, respectively; while the last two terms in Eq. (4) are the familiar Jaynes-Cummings interaction, 
describing the resonant Raman coupling between the two levels ]1) and ]2), which results from the cooperation of the 
cavity mode and the pulse. 

For the case of O = 5, the initial states |2) ]1)^ and ]1) ]0)^ of the whole system, under the Hamiltonian (4), evolve 
as follows 

12)]1), ^ e^«^*/^= [cos(5VAc)|2)|l)c-zsin(5VAc)|l)|0)J, 

|1) |o)^ ^ e^s^*/^= [-isin (5^i/Ac) |2) [1), + cos {gH/\,) |1) 10)J , (5) 

where |0)^ and |1)^ are the vacuum state and the single-photon state of the cavity mode, respectively. The state 
|0) |0)^ remains unchanged under the Hamiltonian (4). 

B. System-cavity-pulse off-resoneint Raman coupling 

Consider n qubit systems labeled by 2, 3, and n + 1. Each qubit system has a four-level structure, as described 
above. The cavity mode is coupled to the |2) o |3) transition of each qubit system, but decoupled (highly detuned) 
from the transition between any other two levels [Fig. 4(b)]. In addition, a classical pulse is applied to each one of 
qubit systems (2, 3, n + \), which is coupled to the jl) -H- |3) transition but decoupled from the transition between 
any other two levels [Fig. 4(b)]. Each pulse has the same frequency oj' . The Hamiltonian for the whole system in the 
Schrodinger picture is 



fe=2 



^ El ]/), {l\ + hg{a+a^,^, + H.c.) + Mlkie'^'^a^,,^ + H.c] 



(6) 



where the subscript k represents the fcth qubit system, Qk is the Rabi frequency of the pulse applied to the kth. qubit 
system, ctjs k ~ \'^)k ' '^13 k — \^)k (^l i and g is the coupling constant between the cavity mode and the |2) o |3) 
transition of each qubit system. 

The detuning between the |1) O |3) transition frequency of the fcth qubit system and the frequency of the pulse 
applied to the kth qubit system is A' = ^31 - co' {k = 2, 3, n + 1), which is identical for each qubit system (i.e., 
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FIG. 4: (a) System-cavity-pulse resonant Raman coupling for a qubit system in a cavity. Here, A = UJ31 — uj is the detuning 
between the pulse frequency uj and the |1) o |3) transition frequency a;3i, while Ac = ujz2 — <^c is the detuning between the 
cavity mode frequency uJc and the |2) -o- |3) transition frequency 0132. The detunings Ac and A are set to be equal, i.e., Ac = A, 
in order to establish the resonant Raman coupling between the two levels |1) and |2) of the qubit system, (b) System-cavity- 
pulse off-resonant Raman coupling for n qubit systems (2, 3, n -|- 1) in a cavity. For simplicity, we here only draw a figure 
for qubit system fe interacting with the cavity mode and a classical microwave pulse (fe = 2, 3, n -\- 1). Here, (5 = Ac — A' 
is the detuning of the cavity mode with the pulse. A' = wai — w' is the detuning between the pulse frequency uj' and the 
|1) -O- |3) transition frequency of qubit system k. Note that 5, Ac, and A' are the same for all qubit systems (2,3, ...,n-\- 1) and 
that the pulses applied to the individual qubit systems (2, 3, n-\-l) have the same frequency uj' . The Rabi frequency of the 
pulse applied to qubit system k is denoted by Q,k. (c) System-cavity off-resonant interaction for n qubit systems (2, 3, n+1) 
in a cavity. The cavity mode is off- resonant with the |2) |3) transition of each qubit system, with a detuning Ac. (d) 
System-pulse resonant interaction. The pulse with a frequency ujr is resonant with the |1) -f^ |2) transition. 

independent of k) because of the same frequency w' of each pulse [Fig. 4(b)]. Under the condition Ac » 5 and 
A' ^ max {122, f^3, ^n+i}, the effective Hamiltonian in the interaction picture can be written as [43] 



"+1 r^2 „2 



^|l)fc(l| + ^«+«|2)fe(2| 



fe=2 

+Xfe(e-^^*a+a+_, + H.c.) 



(7) 



where cr^ ^. = |2)^ (Ij , ^ = Ac — A', and 



Xfe = r!feg(l/Ac + l/A')/2. (8) 
For S :$> g'^/Ac, 6 » max {f2|/A', f2§/A', fi^_,_i/A'}, and S » max{x2,X3) ■■■^Xn+i}, there is no energy exchange 
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between the qubit systems and the cavity mode. Thus, the effective Hamiltonian can be further written as [43-45] 

I — o V c / 



k=2 

n+1 

fc=2 

n+1 



Xk 



(a+a |2),(2|-aa+|l),(l|) 



, XkXk' ( + - - + \ 



(9) 



where the two terms in the second line above describe the photon-number dependent Stark shifts induced by the 
off-resonant Raman couphng, and the two terms in the last parentheses describe the "dipole" coupling between the 
two qubit systems (fc, k') mediated by the cavity mode and the classical pulses. In the case when the level |1) of each 
qubit system is not populated, the Hamiltonian (9) reduces to 



n+l 



k=2 



+ 



xl 



a+a |2),(2|. 



The time-evolution operator for the Hamiltonian (10) is given by 

U{t)=®ltlUkc {t), 



(10) 



(11) 



where Ukc (t) = cxp[i (,g^/Ac + xt/^) \'^)k (^1 1] is the time-evolution operator acting on the cavity mode and the 
/cth qubit system, which takes the following form 



Ukc (t) 



/ 1 \ 

10 

1 

\^ e*^''(*) J 



(12) 



in the basis states 10)^ |0), = (1, 0, 0, 0)^ , 10)^ |1), = (0, 1, 0, 0)^ , |2)^ |0), = (0, 0, 1, 0)^ , and |2)^ |1)^ = (0, 0, 0, 1)^ . 
Here, ipk (t) = exp [i (g^/Ag + xl/^) • Note that the conditional phase (pk is adjustable by varying the Rabi fre- 
quency fifc [see Eq. (8)]. 

C. System-cavity off- resonant interaction 

Consider n qubit systems (2, 3, n -I- 1) with the four-level structure described above. The cavity mode is coupled 
to the 1 2) o 1 3) transition of each qubit system but decoupled (highly detuned) from the transition between any other 
two levels [Fig. 4(c)]. The Hamiltonian for the whole system is given by 



n+l 3 



n-l-1 



H = hw^a+a + E E + ^E '?("^'^23,fc + 



(13) 



fe=2 1=1 



k=2 



where the subscript k represents the fcth qubit system, ctjs k ~ |2)fe ("^1' ^^'^ ^ = 0) 1) 2, 3. 
In the interaction picture, the Hamiltonian (13) becomes 



n+l 

H = fi^5(e-^^=*a+<T23,,+H.c.), (14) 

k=2 

For the case of Ac » 5 (i.e., the cavity mode is off-resonant with the |2) o |3) transition of each qubit system), 
no energy exchange occurs between the qubit systems and the cavity mode. Thus, the effective Hamiltonian can be 
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written as [44,45] 



n+l 2 



fc=2 

n+l 2 

+ ^ J2 (^23,fc^23,fc' +<^23,fc^23,fc') ' 

kjtk'=2 " 

(15) 

where the two terms in the first line above represent the photon-number-dependent Stark shifts, while the two terms 

in the second hne above describe the "dipole" couphng between the two qubit systems (fc, k') mediated by the cavity 
mode. When the level |3) of each qubit system is not excited, the Hamiltonian (15) reduces to 

Hf,fi = -h'^—a^aa22,k- (16) 

The time-evolution operator for the Hamiltonian (16) is given by 

U{t) = ^^llUkc it) , (17) 

where U^c (t) = exp[i (g^/A^) a'^aa22,kt] is the time-evolution operator acting on the cavity mode and the A;th qubit 
system. The operator Ukc (t) is expressed as 



Ukc (t) 



/ 1 \ 

10 

1 

yo e*"^'=(*) J 



(18) 



in the basis states |0)J0)^ = (1,0,0,0)^, |0) Jl)^ = (0, 1, 0, 0)^ , |2) JO)^ = (0, 0, 1, 0)^ , and |2) Jl)^ = (0, 0, 0, 1)^ . 
Here, (j)k (t) = exp (ig'^t/Ac) . 

D. System-pulse resonant interaction 

Consider now a qubit system driven by a classical pulse with frequency LUr and initial phase (p. Moreover, assume 
that the pulse is resonant with the |1) -H- |2) transition, i.e., uir = OJ21, where (jJ2i = {E2 — Ei) /h is the transition 
frequency between the two levels |1) and |2) [Fig. 4(d)]. In this case, the interaction Hamiltonian in the interaction 
picture is given by 

if, = ^(0e^>^|l)(2|+H.c.), (19) 

where is the Rabi frequency of the pulse. From the Hamiltonian (19), it is straightforward to show that a pulse of 
duration t results in the following rotation 

|1) ^ cos^t|l)-ie-''^sin^t|2), 

|2) -ie'-^sin^t |1) + cos^t |2) . (20) 

Above we have introduced four types of interaction of qubit systems with the cavity mode and/or the pulses. The 
results presented above will be employed for the gate implementation discussed in next section. 

III. IMPLEMENTATION OF MULTIQUBIT PHASE GATES 

Let us now consider (n + 1) identical qubit systems in a single-mode cavity or coupled to a resonator, which are 
labelled by 1, 2, and n+l, respectively. The qubit systems {1,2, ...,n + 1) each have the four- level configuration 
depicted in Fig. 3. For each qubit system, the two lowest levels |0) and |1) represent the two logical states of a qubit 
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FIG. 5: (color online) Illustration of qubit systems interacting with the cavity mode and/or the microwave pulses for each step 
of operations during the gate performance. In each figure, blue lines represent the level population of qubit systems before each 
step of operations. Figures from top to bottom correspond to the operations of steps (i)~(vi), respectively. The figures on the 
left side correspond to qubit system 1, while the figures on the right side correspond to qubit system k, with fc = 2, 3, n + 1. 
In each figure, Ac — CJ32 — i^c is the detuning between the cavity mode frequency Uc and the |2) -f-*- |3) transition frequency 
aj32. In (a) and (f), A = ojai — cj is the detuning between the |1) -f-*- |3) transition frequency ojai and the pulse frequency u. In 
(c'). A' — ijJ3i ~ u' is the detuning between the |1) -o- |3) transition frequency ojsi and the frequency uj' of the pulse applied to 
qubit system k, and 5 — Ac — A' is the detuning of the cavity mode with the pulse. 



while the two higher-energy intermediate levels |2) and |3) are used for the gate implementation. In the following, 
qubit system 1 acts as a control while each one of the qubit systems (2, 3, n + 1) plays a target role. 

To realize the proposed gate, let us assume that the level |0) of each qubit system is not affected by the cavity 
mode and/or the pulses during the gate operation; and the cavity mode is off-resonance with the |2) o |3) transition 
of each qubit system with a detuning Ac but decoupled (highly detuned) from the transition between any other two 
levels during the entire operation below. As mentioned above, these conditions can be achieved by prior adjustment 
of the cavity mode frequency [37-41] or the level spacings of the qubit systems before the gate operation [1,2,46], or 
by appropriately choosing the qubit systems (i.e., atoms) to have the desired level strticture. In addition, assume that 
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quantum information is initially stored by the two lowest levels |0) and |1) of each qubit system and the cavity mode 
is initially in the vacuum state |0)^ . 

The procedure for realizing the {n + l)-qubit phase gate described by Eq. (2) is as follows; 

Step (i): Apply a classical pulse to qubit system 1 to induce the Raman transition described in Sec. II A [Fig. 5(a)]. 
It can be seen from Eq. (5) that after a pulse duration tx = ttAc/ (2.g^) , the state |0)^ for qubit system 1 and 
the cavity mode is transformed to the state \2)-^ |l)c • Namely, when the qubit system 1 is initially in the state |1), a 
photon is emitted to the cavity mode after the pulse. On the other hand, the state |0)^ |0)^ remains unchanged during 
the pulse. 

Step (a): Apply a classical pulse (with ujr — 1-021 and (f) — n/2) to qubit system 1 [Fig. 5(b)] and a classical pulse 
(with uij. = LU21 and (j> = — 7r/2) to each of qubit systems (2. 3 n + 1) [Fig. 5(b')]. The duration of each pulse is set 

by t2 = tt/ (2flj (here, fl is the Rabi frequency of each pulse). It can be seen from Eq. (20) that after the pulses, the 

state |2) of qubit system 1 is transformed to the state |1) while the state |1) of each of qubit systems (2, 3, n + 1) 

is transformed to the state |2). Note that since both of the two levels |2) and |3) of qubit system 1 are unpopulated 
after the operation of this step, the qubit system 1 is decoupled from the cavity mode during the operations of steps 
(Hi) and (iv) below [see Figs. 5(c,d)]. 

Step (Hi): Apply a classical pulse (with duration t^) to each of qubit systems (2, 3, r?. + l) to induce the conditional 
phase shift described in Sec. II B [Fig. 5(c')]. Since the level |1) for qubit systems (2, 3, n + 1) is not populated after 
step (ii) above, the effective Hamiltonian of the whole system is given by Eq. (10). Accordingly, the time-evolution 
operator describing this step is the operator given by Eq. (11) for t = ts, i.e., 

U{t3)=(^U2Ukc{t3), (21) 

where Ukc (ts) is the time-evolution operator acting on the cavity mode and qubit system k {k = 2, 3, ...,n + l), which 
takes the matrix of Eq. (12) for t = t^. 

Step (iv): Wait for a time t4. Since no pulse is applied to each qubit system and the cavity mode is off-resonance 
with the |2) ^ |3) transition of each of qubit systems (2,3,...,n + 1) [Fig. 5(d')], this is the case discussed in Sec. II 
C. Note that the level |3) of each qubit system is not excited. Thus, the effective Hamiltonian of the whole system is 
given by Eq. (16). Consequently, the time-evolution operator describing this step is the operator given by Eq. (17) 
for t = t4, i.e., 

C/(i4)=®^±2't^fec(i4), (22) 

where Ukc (^4) is the time-evolution operator acting on the cavity mode and qubit system k, taking the matrix of Eq. 
(18) for t = ti. 

Combining Eq. (22) with Eq. (21) leads to 

|2)fc|l)e ^ exp[ig2(t3+i4)/A,] 

®exp[ix^t3/5] |2)J1),. (23) 

When tz+ti = 2m'Ki^clg'^ {'m is an integer), the transformation (23) becomes 

|2)Jl)e^exp(i^fe)|2)Jl),, (24) 

with 

^'-t*'-^[a; + a^) ^^^^ 

Note that g, 5, Ac, and A' are the same for all qubit systems (2, 3, ...,n + 1). Therefore, the angle 9k can be adjusted 
by changing the Rabi frequency flk of the pulse applied to the qubit system k. The result (24) shows that after the 
operations of steps (i)~(iv), a phase shift exp {i9k) happens to the state |2) of the qubit system k, if and only if the 
cavity mode is in the single-photon state |1)^ . 

Step (v): Apply a classical pulse (with = W21 and </> = — 7r/2) to qubit system 1 [Fig. 5(e)] and a classical pulse 

(with (jjr = 1^21 and = 7r/2) to each of qubit systems (2, 3, n + 1) [Fig. 5(e')]. After a pulse duration t2 = n/ ^2rij , 

the state |1) of qubit system 1 is transformed to the state |2) while the state |2) of qubit systems (2,3, ...,n + 1) is 
transformed back to the original state |1) . 

Step (vi): Repeat the operation of step (i) [Fig. 5(f)]. After a pulse duration ti = ttA^/ (2,g^), the state |2)-^ |1)^ of 
the qubit system 1 and the cavity mode is changed to the state |1)^ |0)^ , which shows that the state |2) of the qubit 
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system 1 is transformed back to the state |1) and the cavity mode returns to its original vacuum state. However, the 
state |0)j^ |0)^ remains unchanged during the pulse. 

One can check tliat the (n + l)-qubit phase gate of one qubit simultaneously controlling n qubits, described by 
Eq. (2), was obtained with (n + 1) qubit systems [i.e., the control qubit system 1 and the target qubit systems 
(2, 3, n + 1)] after the above manipulation. 

To better see how the multiqubit phase gate described by Eq. (2) is implemented by the operations above, let 
us consider a three-qubit example. One can check that the states of the whole system after each step of the above 
operations are summarized below: 



|100)|0), |200)|1), |100)|1), 

|101)|0)^ swi) |201)|1)^ stepfii) |102)|1)^ 

|110)|0), |210)|1), |120)|1), 

|111)|0), |211)|1), |122)|1), 



|100)|1), |200)|1), |100)|0), 

e'^=*|102)|l)^ stepfv) e'^^|201)|l)^ stepfvi) e'^^ |101) |0), 

e''''^|120)|l), ~^ e^^2 |2io) ^ e*^^ |110) |0), 

g«e2e'e3 |i22) |1)^ e^'^^e'^^ \2n) e^^^e*^^ |iii) |o)^ 



(26) 



where \ijk) is an abbreviation of the state \j)2 \k)^ of qubit systems (1, 2, 3) with k S {0, 1, 2}. 
On the other hand, it is obvious that the following states of the whole system 



1000) |0),, 1001) |0),, 1010) |0),, 1011) |0), (27) 

remain unchanged during the entire operation. Hence, it can be concluded from Eq. (26) that a three-qubit phase 
gate of one qubit simultaneously controlling two qubits, described by Eq. (2) (with n ~ 2), was achieved with three 
qubit systems (i.e., the control qubit system 1 and the two target qubit systems 2 and 3) after the above process. 

During the operation of step (ii) or step (v), a single photon is populated in the cavity mode and the state |2) of each 
qubit system is occupied. Therefore, for these two-step operations, there is an accumulated phase shift exp(i2g^t2/Ac) 
to the state |2) of each qubit systc^n. However, when 2t2 <C is + 14, this unwanted phase shift is sufficiently small and 
thus can be neglected. Note that t2 = ir/il and -\-ti = 2mi:/S.c/ 9^ ■ Therefore, the condition 2t2 <C is -|- ^4 turns 
into 17 ^ I (mAc), which can be met by increasing the pulse Rabi frequency ft (i.e., by increasing the intensity of 
the resonant pulses). 

Finally, we should mention that the method proposed here does not work for qubit systems with three levels. 
Suppose; now that the quantum information is encoded in the states |1) and |2) . In this case, during the pulse applied 
to the control qubit system for creating a cavity photon, one cannot have the target qubit systems decoupled from 
the cavity mode. 

IV. TWO TYPES OF MULTIQUBIT PHASE GATES 

Above we have shown how to realize the {n + l)-qubit phase gates described by the operator (2). In this section, 
we focus on two types of multiqubit phase gates, which are depicted in Fig. 6(a) and Fig. 7(a), respectively. It can be 
seen from Fig. 6(a) and Fig. 7(a) that the two types of multiqubit phase gates considered here consist of n two-qubit 
controUed-phase gates each having a shared control qubit but a different target qubit. Below we give a discussion of 
their implementation and importance in quantum information processing. 

A. Implementation of gates of the first-type 

The controlled-phase gate with n-target qubits shown in Fig. 6(a) has the property that the phase for the state |1) 
of each target qubit is shifted by the same amount tt, when the control qubit (i.e., qubit 1) is in the state |1) . Hence, 
this multiqubit phase gate is a special one for Ok = -k above and thus the procedure for realizing it is the same as 
that discussed in the previous section. To obtain this multiqubit phase gate, one needs to set the Rabi frequencies 
of the pulses applied to the target qubit systems (2,3,...,n + 1) to be identical (i.e., Q.2 = ^3 = ■■■ = ^n+i, leading 
to X2 = X3 = ••• = Xn+i = X ) a-iid set the interaction time ts above to be ts = wS/x^, such that 9k = Xk^a/^ = ^ 
(fc = 2,3,...,n + l). 

Note that a controlled-NOT gate with n-target qubits, shown in Fig. 6(b), can also be achieved using the present 
proposal. This is because this gate is equivalent to the n-target-qubit controlled-phase gate discussed here, plus a 
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FIG. 6: (a) Circuit for the first type of controlled-phase gate with qubit 1 simultaneously controlling n target qubits (2, 3, n+ 
1). Here, Z represents a controlled-phase shift e"^ on each target qubit. Namely, if the control qubit 1 is in the state then the 
state |1} at each Z is phase-shifted by tt (i.e., 1 1) — > e^^ 1 1) = — 1 1)), while the state |0) remains unchanged, (b) Relation between 
a n-target-qubit controUed-NOT gate and a n-target-qubit controlled-phase gate. The circuit on the left side is equivalent to 
the circuit on the right side. For the circuit on the left side, the symbol © represents a CNOT gate on each target qubit. If the 
control qubit 1 is in the state |1), then the state at each ® is bit-flipped as |1) — > |0) and |0) — >■ However, when the control 
qubit 1 is in the state |0), the state at each ® remains unchanged. On the other hand, for the circuit on the right side, the 
part enclosed in the dashed-line box represents a n-target-qubit controlled-phase gate. The element containing H corresponds 
to a Hadamard transformation described by |0) (1/a/2) (|0) + |1)) while |1) (1/a/2) (|0) - |1)) . 



single-qubit Hadamard gate acting on each target qubit before and after the n-target-qubit controlled-phase gate 
[Fig. 6(b)] 

The controlled-phase (or controUed-NOT) gates with multiple target qubits shown in Fig. 6(a) [or Fig. 6(b)] are 
important in quantum information processing because they can be applied in entanglement preparation [47], error 
correction [48], quantum algorithms (e.g., the Discrete Cosine Transform [49]), and quantum cloning [50]. 

B. Implementation of gates of the second-type 

The controlled-phase gate with n-target qubits shown in Fig. 7(a) has the property that the phase for the state |1) 
of the target qubit k is shifted by 27r/2'^, when the control qubit (i.e., qubit 1) is in the state |1) . Thus, this multiqubit 
phase gate is a special one for 9k — 27r/2'^ above and therefore can be implemented using the procedure introduced in 
the previous section. To realize it, the Rabi frequencies ^2, ^3, and fln+i for the pulses (respectively applied to 
the target qubit systems 2, 3, and n+l) need to be set non-identical, and satisfy the relation Qk+i/^k = l/v^, 
that is, Xk+i/Xk = l/\/2 (k = 2,3, ...,n); and the interaction time is above needs to be set by = (7r/2) (Vxi)- 
this way, we can obtain 9^ = Xk^a/^ = 2^/2*^ [k = 2, 3, n -t- 1). 

For any two-qubit CP gate described by the transformation |00) |00) , |01) -> |01) , jlO) ^- |10) , and |11) -5> 
e}^ [11) , it is clear that the roles of the two qubits can be interchanged. Namely, the first qubit can be either the 
control qubit or the target qubit, and the same applies to the second qubit. One can see that when the second (first) 
qubit is a control qubit, while the first (second) qubit is a target, the phase of the state |1) of the first (second) qubit 
is shifted by e**^ when the second (first) qubit is in the state |1) , while nothing happens otherwise. Thus, the quantum 
circuit in Fig. 7(a) is equivalent to the quantum circuit shown in Fig. 7(b). It can be seen from Fig. 7(b) that each 
of the n two-qubit CP gates has a shared target qubit (i.e., qubit 1) but a different control qubit (i.e., qubit 2, 3, 
or n-f 1). Note that the n successive two-qubit CP gates shown in Fig. 7(b) are key elements in QFT [51,52]. Hence, 
the importance for this second type of multiqubit controlled-phase gate is obvious. 

The above discussion shows that the two types of multiqubit controlled-phase gates above can be implemented 
by appropriately setting the Rabi frequencies of the pulses applied to the qubit systems. In addition, we should 
point out that by adjusting the Rabi frequencies of the pulses (i.e., changing the intensities of the pulses) applied to 
the target qubit systems, other types of quantum controlled-phase gates with one qubit simultaneously controlling 
multiple target qubits, which may have applications in quantum information processing, can in principle be performed 
using this proposal. 
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FIG. 7: (a) Circuit for the second type of controlled-phase gate with qubit 1 simultaneously controlling n target qubits 
(2, 3, n + 1). Here, the element 27r/2'° represents a controlled-phase shift exp (j27r/2'°), which is performed on the target 
qubit k when the control qubit 1 is in the state |1) (A: = 2,3, ...,«.+ 1). Namely, if the control qubit 1 is in the state 
then a phase shift 27r/2* occurs to the state |1) of the target qubit k, while nothing happens otherwise, (b) Circuit for the n 
successive two-qubit CP gates in a quantum Fourier transform, which is equivalent to the circuit in (a). Here, each two-qubit 
CP gate has a shared target qubit (i.e., qubit 1) but a different control qubit (i.e., qubit 2, 3, or, n + 1). The element 2n/2'' 
represents a controlled-phase shift exp (j27r/2'°), which is performed on the target qubit 1 when the control qubit k is in the 
state |1) (fc = 2, 3, n+l). Namely, if the control qubit k is in the state then a phase shift by 27r/2* occurs to the state 
|1) of the target qubit 1, while nothing happens otherwise. 



V. POSSIBLE EXPERIMENTAL IMPLEMENTATIONS 

In this section, we give a discussion on possible experimental implementations. First, as discussed above, the 
condition 

2*2 < ts + t4 

[i.e., O I (mAc)] needs to be satisfied, which can be achieved by increasing the Rabi frequency tl (i.e., by 
increasing the intensity of the resonant pulses). Second, the total operation time 

T = 2*1 + 2*2 + 13 + *4 = TtAc/ff^ + 2mnAc/g'^ + 2-k/Q. 

should be much shorter than: (i) the energy relaxation time of the level |2) (note that the level |3) is unpopulated 
during the entire operation), and (ii) the lifetime of the cavity mode = Q/2tti/c, where Q is the (loaded) quality 
factor of the cavity. These requirements can in principle be realized, since one can: (i) reduce r by increasing the 
coupling constant g, (ii) increase by employing a high-Q cavity so that the cavity dissipation is negligible during 
the operation, and (iii) choose qubit systems (e.g., atoms) with long spontaneous decay time of the level |2) or design 
qubit systems (e.g., superconducting devices) so that the energy relaxation time 72"^ of the level |2) is sufficiently 
long. 

For the sake of definitiveness, let us consider the experimental possibility of realizing a six-qubit controlled-phase 
gate in a quantum Fourier transform (Fig. 9), using six identical superconducting qubit systems coupled to a resonator 
[Fig. 8(a)]. Each qubit system could be a superconducting chargc-qubit system [Fig. 8(b)], flux-biased phasc-qubit 
system [Fig. 8(c)], or flux-qubit system [Fig. 8(d)]. For the present case, qubit system 1 acts as a target instead 
of a control, while qubit systems (2, 3, 4, 5, 6) play the role of controls (Fig. 9). Each step of operations for 
implementing this gate is the same as that introduced in Sec. III. As discussed above, to implement this gate, the Rabi 
frequencies fl2, ^3, i^i, ^5, and fie of the pulses (respectively applied to qubit systems 2, 3, 4, 5, and 6) need to be set by 
fifc+i/Ofc = 1/V2, leading to Xk+i/Xk = [k = 2, 3, 4, 5). By choosing = lOg, A' = 1002, 02 ~ 0.9^, ~ lOg 

and m = 3, we have: (i) 5 ~ lOg^/Ac lOfi^/A' - 10x2, (ii) 2*2 = tt/O - tt/ (5g) ; and (iii) + t4 - 607r/g, 
where is = (7r/2) {6/X2) ~ 50n/g. Therefore, the condition 2*2 "C *3 + *4 can be well satisfied. Furthermore, for the 
parameters chosen here, the total operation time required for the gate implementation would be r ~ TOw/g. 

As a rough estimate, assume g/27r ~ 220 MHz, which could be reached for a superconducting qubit system coupled 
to a one-dimensional standing- wave CPW (coplanar waveguide) transmission resonator [53]. For the g chosen here. 
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FIG. 8: (color online) (a) Schematic diagram of the setup for six superconducting qubit systems (red dots) and a (grey) 
standing-wave quasi-one-dimensional coplanar waveguide resonator. The two blue curved lines represent the standing wave 
magnetic field, which is in the z-direction. Each qubit system (a red dot) could be a superconducting charge-qubit system shown 
in (b), flux-biased phase-qubit system in (c), and flux-qubit system in (d). The qubit systems are placed at locations where 
the magnetic fields are the same to obtain an identical coupling constant g for each qubit system. The superconducting loop 
of each qubit system, which is a large square for (b) and (d) while a large circle for (c) , is located in the plane of the resonator 
between the two lateral ground planes (i.e., the x-y plane). The external magnetic fiux $ applied to the superconducting loop 
for each qubit system is created by the magnetic field threading the superconducting loop; E,j is the Josephson junction energy 
(0.6 < a < 0.8) and Vg is the gate voltage; A is the wavelength of the resonator mode, and L is the length of the resonator. 
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FIG. 9: Circuit for a six-qubit controUed-phase gate for the quantum Fourier transform. Here, qubits 2, 3, 4, 5, and 6 (on the 
filled circles) play the role of controls acting simultaneously on the target qubit 1. The element 27r/2'^ represents a controUed- 
phase shift exp (i27r/2*) , which is performed on the target qubit 1 when the control qubit k is in the state |1) (fe = 2,3,4,5,6). 
Namely, if the control qubit k is in the state |1), then a phase shift 27r/2'° happens to the state |1) of the target qubit 1, while 
nothing happens otherwise. 

we have r ~ 0.16 /xs, much shorter than '-^ \ [is [2,10]. In addition, consider a resonator with frequency ~ 3 
GHz (e.g., Ref. [24]) and Q ^ 10^, we have ~ 5.3 /xs, which is much longer than the operation time r here. 
Note that superconducting coplanar waveguide resonators with a quality factor Q > 10^ have been experimentally 
demonstrated [54]. 

How well this gate would work needs to be further investigated for each particular experimental set-up or imple- 
mentation. However, we note that this requires a rather lengthy and complex analysis, which is beyond the scope of 
this theoretical work. 

VI. COMPARISON WITH PREVIOUS WORK 

Several points related to this work need to be addressed. First, the present work deals with the realization of 
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multiqubit phase gates for which the phase shifts on each target qubit are tunable, as shown above. Therefore, it is 
much more general than our previous work [55]. Note that in Ref. [55], a method was proposed for the reahzation of 
a multiqubit phase gate which induces a fixed phase-shift of n on each target qubit when the control qubit is in the 
state |1) . The main advantage for the method in [55] is that the cavity mode can be initially in an arbitrary state 
and thus no preparation of the initial state for the cavity mode is required. Second, Ref. [56] showed how multi-ion 
GHZ entangled states can be realized based on performing a multiqubit CNOT gate with one qubit simultaneously 
controlling n target qubits. Since their scheme works for a special case (i.e., a CNOT on each target qubit or a fixed 
phase-shift tt on each target), it fails in performing other types of multiqubit phase gates for which the phase shift on 
each target qubit is not tt. Third, we note that using the method presented in [56], a multiqubit CNOT gate or phase 
gate with a phase shift n on each target qubit (i.e., the first type of multiqubit gate discussed above) is rather difficult 
to implement in the other systems due to different physical mechanisms. Last, the present proposal for realizing the n 
successive two-qubit CP gates in QFT is different from the previous work reported in [57]. As shown above, using the 
present proposal, the n successive two-qubit CP gates in QFT can be simultaneously performed using one single-mode 
cavity only. However, as argued in Ref. [57], when using the proposal in [57] to implement n successive two-qubit 
CP gates in QFT, n single-mode cavities would be required or a single cavity with various modes interacting qubit 
systems (e.g., atoms) needs to be designed. 

VII. CONCLUSION 

We have presented a method to realize a multiqubit tunable phase gate with one qubit simultaneously controlling 
n target qubits in a cavity QED. As shown above, the method has the following features: (i) Neither adjustment of 
the level spacings of qubit systems nor adjustment of the cavity mode frequency during the gate operation is needed, 
thus the operation is much simplified; (ii) The operation time required for the gate realization is independent of the 
number of qubits and thus does not increase with the number of qubits; (iii) The n two-qubit CP gates involved can 
be simultaneously performed, which significantly reduces the gate operation especially when the number n of qubits 
is large; (iv) The phase shift on the state of each target qubit can be adjusted by changing the Rabi frequencies of the 
pulses applied to the target qubit systems; and (v) Certain types of significant quantum controUed-phase gates with 
one qubit simultaneously controlling multiple target qubits or quantum controUed-phase gates with multiple control 
qubits simultaneously acting on one target qubit (e.g., the multiqubit phase gate in the QFT discussed above) can be 
performed by using this proposal. 
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APPENDIX: HOW TO HAVE THE LEVEL |0) NOT TO BE AFFECTED DURING THE 

OPERATION 

As shown above, three types of interaction are needed for the gate implementation, which are the system-cavity- 
pulse resonant Raman coupling, the system-cavity-pulse off-resonant Raman coupling, and the system-pulse resonant 
interaction. For the last type of interaction, since a resonant pulse is applied, the level |0) can be easily made not 
to be affected by the pulse, by prior adjustment of level spacings such that the transition between the level |0) and 
any one of other levels is largely detuned (decoupled) from the pulse. In the following, we will focus on how to have 
the level |0) not to be affected by the pulse and the cavity mode, for the case when the pulse is off-resonant with the 
|1) -H- |3) transition and the cavity mode is off-resonant with the |2) -o- |3) transition, i.e., the case which applies to 
the former two types of interaction. 

To begin with, we define Ai = \{Ei - Eq) /h - a;| , A2 = [(£'2 - £^0) - ^1 , and A3 = [(£3 - Eq) /h - uj\ as the 
large detuning of the pulse with the |0) •<-> |1) transition, the large detuning of the pulse with the |0) -O- |2) transition, 
and the large detuning between the pulse and the |0) o |3) transition, respectively (Fig. 10). Here, oj is the pulse 
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FIG. 10: (Color online) Large detuning of the pulse/the cavity mode with the transition between the level |0) and any one of 
other levels. 



frequency. In addition, we define A J — \ (E\ ~ E^) jfi — liJc\ and ~ |(£^2 — £^i) /?i — wd as the large detuning of the 
cavity mode with the |0) |1) transition and the large detuning of the cavity mode with the |1) ^ |2) transition, 
respectively (Fig. 10). 

A. Charge-qubit system 

A charge-qubit system has the level structure depicted in Fig. 3(a) [or Fig. 10(a)]. From Fig. 10(a), one can see that 
the dotted line falls within the range between the levels |1) and |2) in the case when the pulse is off-resonant with the 
|1) O 1 3) transition. This is because the level spacing E2 — Eq is larger than the level spacing E^ — Ei, resulting from 
E3 — E2 < El — Eg as mentioned earlier [see the caption of Fig. 3(a)]. Fig. 10(a) also shows that the dot-dashed line 
falls within the range between the levels |0) and 1 1) in the case when the cavity mode is off resonant with the 1 1) O |3) 
transition. This is due to E3 — E2 < Ei — Eq. In order to have the level |0) not to be affected by the pulse and the 
cavity mode, one will need to adjust the level spacings of the qubit system before the gate operation to achieve: (i) 
the large detunings Ai and A2 [Fig. 10(a)], such that neither |0) O |1) transition nor |0) o |2) transition is induced 
by the pulse; and (ii) the large detuning [Fig. 10(a)] such that the cavity mode does not cause the |0) o |1) 
transition. From Fig. 10(a), it can be seen that the large detuning regime for the pulse and the |0) O |3) transition 
is automatically satisfied, and the large detuning of the cavity mode with the |0) ^ |2) or |0) O |3) transition is also 
automatically met. Therefore, by prior adjustment of the level spacings to have the large detunings Ai, A2, and Aj, 
the level |0) will not be affected by the pulse and the cavity mode during the operation. 

B. Phase-qubit system 

A phase-qubit system has the level structure shown in Fig. 3(b) [or Fig. 10(b)]. It can be seen from Fig. 10(b) 
that when the pulse is off resonant with the |1) |3) transition, the dotted line falls within the range between 
the two lowest levels |0) and |1), which could be achieved by adjusting the level spacings of the qubit system. In 
addition. Fig. 10(b) shows that the dot-dashed line falls within the range between the levels |1) and |2), in the case 
when the cavity mode is off resonant with the transition between the upper two levels |2) and |3). This is because of 
E2 — El > £'3 — E2 [see the caption of Fig. 3(b)]. To have the level |0) not to be affected by the pulse and the cavity 
mode, one will need to adjust the level spacings of the qubit system before the gate operation to obtain: (i) the large 
detuning Ai [Fig. 10(b)] such that the |0) O |1) transition induced by the pulse is avoided; and (ii) the large detuning 
A^ which ensures the large detuning between the cavity mode and the |0) ^ |1) transition due to Ei — Eq > E2 — Ei 
[Fig. 10(b)]. In addition. Fig. 10(b) shows that the large detuning regime for the pulse and the |0) |2) or |0) O |3) 
transition is automatically met, and the large detuning of the cavity mode with the |0) o |2) or |0) o |3) transition 
is automatically satisfied. Thus, as long as the large detunings Ai and A^ are met by prior adjustment of the level 
spacings, the level |0) will not be affected by the pulse and the cavity mode during the operation. 
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C. Flux-qubit system 

A flux-qubit system has the level structure in Fig. 3(c) [or Fig. 10(c)]. Fig. 10(c) shows that the dotted line falls 
within the range between the levels |2) and |3) in the case when the pulse is off-resonant with the |1) ^ |3) transition. 
This is because the level spacing E2 — Eq is smaller than the level spacing E3 — E\ due to Ei — Eq < E^ — E2 [see 
the caption of Fig. 3(c)]. Fig. 10(c) also shows that the dot-dashed line is above the level |1) , which could be reached 
by adjusting the level spacings of the qubit system. In order to have the level |0) not to be affected by the pulse and 
the cavity mode, one will need to adjust the level spacings of the qubit system to achieve: (i) the large detunings 
A2 and A3 [Fig. 10(c)], such that neither |0) ^ |2) transition nor |0) O |3) transition is induced by the pulse; and 
(ii) the large detuning [Fig. 10(c)], in order to avoid the |0) O |1) transition induced by the cavity mode. From 
Fig. 10(c), it can be seen that the large detuning regime for the pulse and the |0) ^ |1) transition is automatically 
satisfied, and the large detuning of the cavity mode with the |0) ^ |2) or |0) |3) transition is automatically met. 
Hence, as long as the large detunings A2, A3, and Aj are satisfied by prior adjustment of the level spacings, the level 
|0) will not be affected during the operation. We here should mention that because the dotted line falls within the 
range between the upper two levels |2) and |3) [Fig. 10(c)], a more careful adjustment of the level spacings would 
be needed to obtain the large detunings Ai and A2 for the flux-qubit system, when compared with the charge-qubit 
system or the phase-qubit system. 
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